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We propose a systematic way to investigate the low- 
temperature thermodynamic properties of quantum spin sys- 
tems subject to the restriction that only a finite number of 
bosons may occupy a single lattice site. Such a kinematical 
interaction results in appearance of a temperature dependent 
chemical potential. Its low-temperature asymptotics is calcu- 
lated self-consistently using the functional integration tech- 
nique. 
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To systematically calculate the thermodynamic prop- 
erties of a two-dimensional (2D) quantum ferromagnet 
at low temperatures remains an unsolved problem of the 
spin wave theory The main difficulty is in how to 

consistently calculate both dynamical and kinematical |^] 
interactions of magnons. Such interactions are believed 
to vanish the average spin (magnetization) of a system, 
as it naturally follows from the absence of long range 
order in 2D systems at low temperatures. 

Recently, a few different approaches to this problem 
have been developed. In his "modified spin wave the- 
ory" , Takahashi Q uses the Holstein-Primakoff represen- 
tation to introduce (by hands) a chemical potential for a 
system of bosons. This chemical potential is essentially 
nonzero in 2D systems, even in the absence of magnetic 
field. Thus, the zero magnetization condition, {Sz) — 0, 
is enforced. Despite good agreement with numerical cal- 
culations, such an approach seems to be not ultimately 
self-consistent. In another approach by Arovas and Auer- 
bach a functional integral for a partition function is 
constructed using the generalized bosonic SU (N) repre- 
sentation of spin algebra. Being calculated using the 1 /N 
expansion around the mean field {N — oo) saddle point, 
the results were then extended to the case of = 2 [i.e., 
SU{2)]. The validity of this procedure has not, however, 
been justified. 

In this letter, we propose a self-consistent method to 
calculate the partition function of a system of bosons 
with kinematical interaction, i.e., the number of bosons 
on a lattice should not exceed some number L. One 
should note that in contrast to Takahashi Q, who con- 
structed a theory of low-dimensional ferromagnet by 
means of introducing "by hand" a chemical potential 
for Holstein-Primakoff bosons to satisfy the condition 
(Sz) = 0, we develop a theory with the chemical potential 
as a result of controllable approximations. We write the 
partition function in the functional integral representa- 
tion. In the one-loop approximation, we then show that 
the system of kinematically interacting magnons is equiv- 
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FIG. 1. A lattice site; up to L bosons may occupy this site. 



alent to that of non-interacting magnons with some non- 
zero, temperature dependent chemical potential. Its low- 
temperature behavior is finally calculated. One should 
comment that the theory of the systems with kinematical 
interaction which limits the maximum number of parti- 
cles per site is, in fact, closely related to the parastatistics 
theory ^ with essentially similar "exclusion" principles. 

We demonstrate the method on a simple model with 
the Hamiltonian: 



H = ^t^k^k^k , t^k 
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on a square lattice, where a is a lattice constant, (3—1/1 
is the inverse temperature, and 6 and 6+ are the Bose 
operators. 

As mentioned above, we consider the system of bosons 
with kinematical interaction. In other words, an each 
lattice site can be occupied by no more than L bosons. 
Thus, we need a projecting operator to eliminate all un- 
physical states with n > L + 1 (see Fig. I), while all 
physical states < n < L are to remain unchanged. 
Heuristically, we define the projecting operator P — 
Y\ - Pi as follows: 



P 



L 



(2) 



where : O : is normally ordered form of an operator 0. 
Here we used a useful identity: 



which can be proved straightforwardly. 

Now the partition function of the system is: 

Z = Sp{Pe-f^"}. 
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To obtain the functional representation for the partition 
function (Q) , it is instructive to transform to the coherent 
basis 0: 



\u) = e "2"^ 

with the properties: 

J\u){u\^-^^Y.n=,W){n\^i 
b\u) = u\u) , — {u\u* , 

(v\u) = e"^("'"+''''')+''*" . 



(5) 
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Using identity (pa]), we rewrite the partition function as 
follows: 



[v\P\u){u\e-'^"\v)DuDv. 



(7) 



Since all the operators are in normal ordered form, we 
use Eqs. (6b) to obtain the following matrix elements: 



n=0 



{u\e-P"\v) = e 



Thus, the exact partition function becomes: 

Z = j e-'^^'^'^'^DuDv, 
where the "action" S{u,v) is 

S{u, v) = ^(ufcUk + Wk^k - e~''"''WkUk) 
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In the mean field approximation, we expand the loga- 
rithm up to the first nonlinear term: 
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(11) 



It can be shown that the higher order terms omitted in 
Eq. ( p^ ) give only small corrections to the expressions 
obtained at low temperatures, because of the small pa- 
rameter |logT|~^ <C 1. Now the action S{u,v) reads 
as 



S{u, v) « Soiu, v) + Sintiu, v) 



(12) 



where 



Sq{u, v) = ^(uk"k + t'k^'k - e ^"''UkWk - Uk«k) , (13a) 

(13b) 



The term Sint{u,v) is of fundamental importance. In- 
deed, omitting it from Eq. ( p^ yields the partition hmc- 
tion: 



Zo = / e-^°("'")DuDw, 



which suffers from infrared divergences in the thermody- 
namic average for the one- and two-dimensional cases: 



00 



Thus, it is equivalent to the omission of the projection 
operator from Eq. (^. 

We now linearize the "action" Eq. (|l|) as follows: 

S^nt{u, V)^Y1 ((^ + 1)^''<"' - LA^ + ') (14) 



and introduce the mean field partition function 

Zmf{A)= J e-^^f^"^'''^'>DuDv , (15) 

with 

Smfiu, V, A) = ^(u^Uk + v^Vk - e^'^'^'^v^uu 

k 

- (1 - (i + l)A^)u;;?;k - iA^+i) . (16) 

Here we defined the thermodynamic average: 

A = {v*u,)^f . (17) 

Thus, in the mean field approximation, we have the 
quadratic mean field "action" Smfiu^v, A) instead of 
nonlinear the "action" S{u,v). The parameter A will 
be self-consistently obtained below, in accordance with 
Eq. (0). The simphcity of the mean field "action" (|l|) 
allows us to perform functional integration in the parti- 
tion function to yield: 



-^-'^(l- (i + l)AL) 



(18) 



Here 

F,„,(A) = i5:(-LA^+i 



+ In (1 - e-'^"''(l - (L + 1)A^))) . (19) 
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(25) 



0.05 



L»l 




0.00 



FIG. 2. The temperature dependence of Ao for two cases: 
L = 1 and L > 1. 



From the condition: 

dFrafiA) 



dA 



= 0, 



we obtain the self-consistent equation for Aq 

1 



Ao = -V 



N ^ - (1 - (i + 1)A^; 



(20) 



(21) 



In the two-dimensional case at low temperature (and for 
Ao ^ 1), it reduces to 



Ao = 



T ln((£ + l)Ao^) 
2ttJ 1 - (L-h 1)A^ 



-^ln((X + l)A^) 



(22) 



The temperature dependence for Aq is presented in Fig. 
H for two limiting cases: L — 1 and L ^ 1 (here L — 
500). The thermodynamic average Aq vanishes at T = 
and rapidly grows as temperature increases. When = 
r/27rJ< 1, 



An 



-T,i In (T*) 



(23) 



Now, it is obvious that we can exclude the project- 
ing operator from Eq. (^), but simultaneously we must 
introduce the chemical potential 



= T{L + 1)A^ = 2^g-27rJAo/T 



(24) 



in Eq. (|l|). This is precisely the Takahashi's assump- 
tion that results in very good agreement with the Bethe- 
anzatz for an one-dimensional ferromagnet at low tem- 
perature. Note, n decreases with L and vanishes in the 
absence of the kinematical interaction, L = oo (i.e., non- 
interacting bosons). From Eqs. (23), (24), we obtain the 
low-temperature asymptotics: 



To conclude, we have shown that a kinematical inter- 
action in 2D magnon systems can be compensated via 
introducing a "fictitious" chemical potential into their 
effective dispersion law, and proposed a systematic, self- 
consistent procedure of calculating its low-temperature 
dependence in the one-loop approximation. 

We would like to acknowledge P.H. Diamond for his 
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